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Abstract. We study dyonic soliton and black hole solutions of the su(2) Einstein- 
Yang-Mills equations in asymptotically anti-de Sitter space. We prove the existence 
of non-trivial dyonic soliton and black hole solutions in a neighbourhood of the trivial 
solution. For these solutions the magnetic gauge field function has no zeros and we 
conjecture that at least some of these non-trivial solutions will be stable. The global 
existence proof uses local existence results and a non-linear perturbation argument 
based on the (Banach space) implicit function theorem. 



1. Introduction 

The study of soliton and black hole solutions of Einstein- Yang-Mills (EYM) theory 
has a long history, sparked by the discovery of four dimensional, asymptotically flat, 
spherically symmetric, solitons [1] and black holes [2] when the gauge group is su(2) and 
the gauge field is purely magnetic. Although these solutions are unstable [3], solutions 
of EYM in asymptotically flat space have been widely studied (see, for example, jl] for 
a review). 

Changing the space-time asymptotics dramatically changes the properties of EYM 
solutions. In this paper we focus on four-dimensional, asymptotically anti-de Sitter 
(adS) space-times. The first EYM solutions in adS were found for gauge group su(2) 
and, as in asymptotically flat space, spherically symmetric black hole [5] and soliton [6] 
solutions exist for a purely magnetic gauge field. However, unlike the asymptotically 
flat solutions, at least some of these purely magnetic, asymptotically adS, solutions are 
stable under linear perturbations [5j El [7J. Similar results hold if the gauge group is 
enlarged to su(iV) [HIE], yielding black holes with abundant gauge field hair [TO] . 

As well as the existence of stable solitons and black holes with a purely magnetic 
gauge field, su(2) EYM solutions in adS have another surprising feature: the existence 
(in numerical simulations) of spherically symmetric dyonic solutions for which the gauge 
field has non-trivial electric and magnetic parts [6]. This is in contrast to the situation 
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in asymptotically flat space-time, where the only EYM solution where the gauge field 
has a non-zero electric part is the trivial embedded Reissner- Nordstrom solution [TTj . 

In this paper we prove analytically the existence of four dimensional, spherically 
symmetric, dyonic soliton and black hole solutions of su(2) EYM in adS, as previously 
found numerically [6]. In section |2] we briefly summarize our model, the field equations 
and boundary conditions. We also review the key features of the numerical solutions [6] . 
We begin our analytic work in section [31 by proving the local existence of solutions of the 
field equations near the origin, event horizon and infinity. In section HJ we then prove 
the existence of global solutions of the field equations, using a non-linear perturbation 
argument based on the implicit function theorem rather than the usual "shooting" type 
methods which were employed in, for example, [5]. Our conclusions are presented in 
section El 



2. su(2) Einstein- Yang-Mills theory 

In this section we present the field equations for su(2) Einstein- Yang-Mills theory in 
adS, and briefly review the properties of the numerical dyonic solutions [6]. 

2.1. Ansatz and field equations 

We consider four-dimensional EYM theory with a negative cosmological constant A < 0, 
with the following action: 

1 



Seym = -J #xyf^ [R-2A- Tr F, W F^\ , (1) 

where R is the Ricci scalar, A the cosmological constant, the Yang-Mills gauge field 
and Tr denotes a Lie algebra trace. We have fixed the gauge coupling constant g = 1. 
Throughout this paper, the metric has signature (— , +, +, +) and we use units in which 
47rG = c = 1. 

The field equations derived from varying the action flTJ are: 

= D^FS = V M F/ + [Ap, F^] , (2) 
where the YM stress-energy tensor is 

Tfu, = Tr FpxF v x - -g^ Tr F Xa F X(J , (3) 

and F^y is given in terms of the gauge potential by 

F^ = d^A v - d u A^ + [Ap, A v ] . (4) 

In 01]), square brackets denote the Lie algebra commutator. 

Considering static, spherically symmetric space-times, the metric in standard 
Schwarzschild-like co-ordinates takes the form 

ds 2 = -/i(r)5(r) 2 dt 2 + Mr)" 1 dr 2 + r 2 d6 2 + r 2 sin 2 9 dip 2 , (5) 
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where the metric functions /i(r) and S(r) depend on the radial co-ordinate r only. The 
metric function /i(r) can be written as 

, , 2m(r) r 2 , . 

tir) = l-—?- L + p, (6) 

where the adS radius of curvature £ is defined by 

'-4 

With a suitable choice of gauge, the su(2) gauge field takes the form [12] 




(7) 



A = Adt+i(C -C H ) ol6- 1 - [{C + C H )sm6 + DcosO] dip, (8) 
where A, C and D are 2x2 matrices, given by 

A -if a(r) ^ r - f \ n - 
2\ -a(r) )> \0 J ' 

with a(r) and w(r) being real functions of r only. We emphasize that our gauge field 
ansatz (jSJ) has both an electric part (Adt) described by the single function a(r) and 
a magnetic part described by the single function u(r). We will refer to a(r) as the 
"electric gauge field function" and u(r) as the "magnetic gauge field function". 

Substituting the metric (jHJ) and gauge field (jBJ) into the field equations (j2J) we obtain 
two Einstein equations: 

, rV 2 aV (u 2 -lf 
y_2«V 2a/ 2 

5 ~ r/i 2 S 2 + r ' ( j 

and two Yang-Mills equations: 



2a' a'S' 2au 2 
a" = + — - + -J—, 

„ uj'S 1 u'fi' a 2 u u(u 2 — l) 

u — jr~^ + • 

In the above equations, a prime ' denotes d/dr. The field equations f TTOMTTT) are invariant 
under the transformations a — > —a and uj — > —uj separately. They are also invariant 
under the following scaling symmetry: 

S(r) -> XS(r), a{r) -)• Aa(r), (12) 

for any constant A. This corresponds to the invariance of the form <^ of the metric and 
the form (jSJ) of the gauge potential under scalings t — >■ \~H. 

Finally we note that, unlike the situation for purely magnetic su(2) EYM [5l [6], 
here the second Einstein equation for the metric function S(r) does not decouple from 
the other equations. 
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2.2. Boundary conditions 

The field equations ( fTUHTT]) have three singular points of interest^, corresponding to the 
origin r = 0, black hole event horizon r = r^ (defined by zeroes of fi, if there are any) 
and infinity r — > oo. In this section we state the boundary conditions on the functions 
a, u, fj, and S at each of these singular points, before proving in section |3] that solutions 
of the field equations satisfying these boundary conditions exist. 

2.2.1. Origin We start by assuming that the field variables have regular Taylor series 
expansions about r = 0: 

m(r) = m + m^r + m 2 r 2 + m 3 r 3 + 0{r ), 

Sir) = S + S ir + S 2 r 2 + 0(r 3 ), 

a(r) = a + cmr + a 2 r 2 + a 3 r 3 + 0(r 4 ), 

u(r) = ujQ + LOir + u 2 r 2 + 0{r z ). (13) 



To avoid a singularity either in the metric (J5J) or the field equations fUOHTT]) . it must be 
the case that 

m = mi = i7i2 = Si = ao = a 2 = oui = (14) 

and u = ±1. We take u = 1 without loss of generality since the field equations are 
invariant under the map u — > —u. The constants 7713, S 2 and a 3 are then determined 
in terms of So, «i and u 2 by the field equations and the series ( ITBl become [6]: 



m(r)= (^ + 2^ r 3 + 0(r 4 



S(r) = So + (j| + 4Sow2) ^ + 0(r 3 ), 

tt (r) = ttl r + ^ ^ + 8a; 2 2 + 2u 2 - r 3 + 0(r 4 ), 

w(r) = l + w 2 r 2 + 0(r 3 ). (15) 

Subject to proving convergence, these series give us a three-parameter family of solutions 
near the origin, parameterized by So 7^ 0, a.\ and u 2 (plus the adS radius of curvature 
£). In practise the value of Sq is fixed by the requirement that S(r) — > 1 as r — > 00. 



2.2.2. Event horizon For a regular, non-extremal black hole event horizon at r = r^, 

the value of m(rh) is fixed by the condition /i(r^) = to be 

3 

= y + ^. (16) 

I There may in addition be singular points corresponding to zeroes of 5". However the second Einstein 
equation (fTU)l is a first order, linear inhomogeneous equation in S 2 showing that, when solutions of the 
full system exist, (S 2 )' > 0. We will be interested in solutions with a positive initial value for S, and 
so S — cannot arise. This is proven more carefully below. 
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The first Yang-Mills equation (ITT]) has a singularity at r = r h unless the electric gauge 
field function a vanishes there. Series expansions of the field functions near the event 
horizon therefore take the form [6]: 

m ( r ) = -J + 2^2 + m 'h ( r - Th ) + °( r ~~ r ^) 2 ' 
S(r) = S h + S' h (r-r h ) + 0(r-r h ) 2 , 
a{r) = a' h (r - r h ) + 0{r - r h ) 2 , 

u(r) = uj h + u' h (r - r h ) + 0{r - r h ) 2 , (17) 
where m' h , S' h and oj' h are given by the field equations in terms of a' h and ojh as: 
rla'l , H - if 



2S2 2r^ 



^ = 2a' 2 ul + 2<g fe 
ft S h r h fi'(r h ) 2 r h 



J °°h ( 2 



M-i), (18) 



and 



„'(r») = I - M + ^ > o. (19) 
r h r h t l 

In this case we also anticipate the existence of a four-parameter family of solutions, 
with parameters r^, Sh, ot'h and Uh (as well as the adS radius of curvature £). As with 
the expansions near the origin, the value of Sh will be fixed by the requirement that 
S(r) — > 1 as r — > oo. 

2.2.3. Infinity As r — > oo, we require our metric (EJ) to approach adS, so that S'(r) — > 1 
and m(r) — > M where M is a constant which corresponds to the mass of the soliton 
or black hole. Assuming that the field variables have regular Taylor series expansions 
near infinity, the second of the Einstein equations ( ITOl) implies that S'(r) = 0(r~ 5 ) as 
r — > oo, so the series take the form 

m{r) =M-±(| + alulf + | + | + 0(r" 

(«L^ 4 + c?)+0(r- 5 ), 



5(r) 






= 1 - 


2r 4 ' 


a(r) 








+ — 




r 


w(r) 


= Woo 


r 



9(r~ 2 ). (20) 

The field equations place no restrictions on the constants a^, di, Wqo, C\ or M, giving 
an expected five-parameter family of solutions (which depend on the adS radius of 
curvature £ as well as the parameters above). 
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2.3. Trivial solutions 

The field equations (fTUI - fTTT) have a few simple solutions: 

(i) If we set 

a(r)=0, w(r) = ±l, (21) 

then both m(r) and S'(r) are constants, giving the Schwarzschild-adS black hole 
solution. If, in addition, m(r) = then pure adS is also a solution of the system. 

(ii) Alternatively, if we set 

a{r) = 0, w(r) = 0, (22) 
then S(r) is still constant, but in this case 

m(r) — M , (23) 

2r 

giving the metric for a magnetically charged, Abelian, Reissner-Nordstrom-adS 
black hole with charge Qm = 1- 

(iii) The Abelian Reissner-Nordstrom-adS black hole with electric charge Qe and 
magnetic charge equal to unity is also a solution, obtained by setting 

a(r) = 2* w = 0, (24) 
r 

in which case 5(r) = 1 and 

m(r) = M- ( Q2js + X \ (25) 

Due to the coupling between the electric and magnetic gauge field functions in the 
Yang-Mills equations ( ITTj) . Abelian Reissner-Nordstrom-adS with non-zero electric 
charge but zero magnetic charge is not a solution of the system. 

(iv) Finally, by setting a(r) = 0, the field equations ( TlTil - fTTI) reduce to those for purely 
magnetic su(2) EYM theory with a negative cosmological constant, whose solutions 
have already been studied in depth [3 El |9]. In this case the second Einstein 
equation for the metric function S(r) (fIU|) decouples from the other equations. 

Our focus in this paper is non-trivial solutions (that is, solutions which do not appear in 
the list above), in a neighbourhood of either the pure adS solution or Schwarzschild-adS 
solution. 



2.4- Properties of the numerical solutions 

The dyonic field equations ( fTUMTTT) were solved numerically in [6], so here we briefly 
review those results. The field equations are solved by starting the integration close 
to either the origin (for soliton solutions) or event horizon (for black hole solutions), 
using the relevant series (( I15j) or (117j) respectively) and integrating outwards to large 
r, stopping either once the field variables have converged to their asymptotic values 
within a desired numerical tolerance or the solutions have become singular. To satisfy 
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Figure 1. Typical dyonic soliton solution with t = 1, a\ = 0.2455 and ui-2 = 0.2. Both 
the electric gauge field function a(r) and the magnetic gauge field function uj(r) are 
monotonically increasing and neither has any zeros for r > 0. 



the boundary condition on S(r) at infinity, it is simplest to initially set either Sq = 1 (for 
soliton solutions) or Sh = 1 (for black hole solutions) and integrate the field equations. 
The value of S(r) as r — > oo will then be Sex, which is in general not equal to unity. 
Performing a scaling transformation (I12p with A = S 1 ^ 1 then gives a numerical solution 
for which the boundary conditions fl20l) hold. With this approach, the parameters 
which can be varied are ot\ and U2 for soliton solutions and r^, a! h and Uh for black hole 
solutions, as well as the adS radius of curvature I. 

For each value of I and (with = for solitons), regular soliton and black 
hole solutions are found in continous regions of the two-dimensional parameter space 
(a.i,uj2) or (a' h ,Uh)- The size of these regions increase as the adS radius £ decreases (or, 
alternatively, as |A| increases). Some typical solutions are shown in figures [1] and [2j 

It is straightforward to show from the field equations (fl^HTTl) that the electric 
gauge field function a(r) is monotonic, as are the metric functions m(r) and S(r). This 
is not the case for the magnetic gauge field function u(r), which typically has some 
zeros before approaching its asymptotic values. However, we find solutions in which 
u(r) has no zeros (that is, is nodeless) in a neighbourhood of the trivial pure adS (for 
solitons) or Schwarzschild-adS (for black holes) solution. We observe that the size of 
the neighbourhood in which these nodeless solutions exist shrinks as the adS radius of 
curvature I increases (or |A| decreases). Similar behaviour is observed in the purely 
magnetic case [SI IB]- Further details of the space of solutions can be found in [Bl [13] . 

3. Local existence 

In this section we prove the existence of local families of solutions of the field equations 
fTTUlTilT) in neighbourhoods of the singular points r = (origin), r = (black hole 
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Figure 2. Typical dyonic black hole solutions with I = 1, ru = 1, oi h — 0.2971 and 
ujh = 0.8. The electric gauge field function a(r) is monotonically increasing and the 
magnetic gauge field function ui(r) is monotonically decreasing. Neither gauge function 
has any zeros for r > r^. 



event horizon) and r — > oo (infinity). We follow the method of [Hj, making use of the 
following theorem: 

Theorem 1 [Lffl Consider a system of differential equations for n + m functions 
u = (ui, u 2 , ■ ■ ■ u n ) T and v = (vi, v 2 , ■ ■ ■ v m ) T of the form 

X ~dx = X ( x > u > v ) ' 

X ~d~x = ~ T ' iVi + X ^ 9i ^' U ' ^ ' ^ 
with constants n > and integers > 1, and let C be an open subset of W n such 

that the functions fi, % = 1, . . . n and g\, % — 1, . . . m are analytic in a neighbourhood of 
x = 0, u = c, v = for all c G C. Then there exists an n-parameter family of solutions 
of the system If2h]) such that 

Ui {x) = c i + 0(x tTi ), Vi(x) = 0(x^), (27) 

where Ui(x), i = 1, . . .n and Vi(x), i = 1, . . .m are defined for c G C, for \x\ < xq(c) 
(for some Xq(c) > 0) and are analytic in x and c. 

The key aspect of applying this theorem is setting the field equations in the 
required form (1261) near each of the singular points. Proving the local existence of 
families of solutions with the desired boundary conditions (jl5| \T7\ 120]) then follows in a 
straightforward way from the above theorem. One advantage of using Theorem [1] is that 
the local solutions are analytic in the independent variable x and also in the parameters 
c. This is helpful when proving the global existence of solutions in section HJ In each case 
the proof of local existence is very similar to that in [H], although the presence of the 
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electric gauge field function a(r) complicates the analysis both in terms of the boundary 
conditions at the origin and event horizon, and also because, unlike the situation for 
purely magnetic configurations, here we must include the field equation for the metric 
function S(r). For each of the three propositions in this section, we also add the equation 
d£/dx = (for each particular choice of independent variable x), which clearly has the 
trivial solution £ = £ . Adding this equation to our system means that the local solutions 
whose existence we prove will also be analytic in £ , or, equivalently, in i. 

Proposition 2 (Local existence of solutions in a neighbourhood of the origin) 

There exists a three-parameter family of local solutions of the field equations (Q2HI2]] 

near r = 0, satisfying the boundary conditions ( f73j) and analytic in S , a\, u 2 , r and £. 

Proof Following [14], we take our independent variable x — r, and we define the 
following dependent variables: 

Ai(x) = ^[l-/i(r)], \ 2 = S(r), 
M*) = —> M*) = \Hr)-l], 

t . . 1 da a t . , adw , . 

6W = ; Tr~7v Ux) = ~rTr (28) 



These dependent variables are the same as in the corresponding proof in [14], apart 
from ipi and £i which depend on the electric gauge field function a(r). In terms of these 
variables, the field equations (JTDHTTJ) take the form 



dXi 3 M\ 

~ " 3Al ~ e + "aT 



X ~JT = " 3A i " 77 + "Jr + 4 ^ I + 2 & + a# 0) i, 



^2 2p 

dx 

x—— = xfi, 
dx 

dtb 2 9 

= -2^ 2 + 6 + X 2 Jq,2, 

dx 

x^p- = -4& + xH ,i, 
dx 

= + 2?p 2 + x 2 n , 2 , (29) 

dx 

where Qo,i and Ho,i are polynomials in the Aj, (1 — a; 2 Ai) \ A 2 1 , the ^ and the and 
are given explicitly by: 

0o,i = § (2^i + z£i) + ^ (4 + x 2 ^) + 



x 2 Ai 



2x^1 (Ai + 2^ 2 + x 2 jj 2 ) 
+ (l-o^A^A 2 
2^(1 + a 2 ^) 2 2A 2 £ 2 
^ (l-x 2 A0 2 A 2 + (l-x 2 A0 2 ' 
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A further transformation is needed to cast the equations ([29]) into the required form 
(126]) . In particular, we define new variables uq, vq and Ai as follows [1] 

u = - {1P2 + 6) , = ^ (2^2 - 6) > 



Ai = A x + - - ^ - 4^ + 2< (31) 



£2 A^ 

in terms of which the equations for Ai, ip 2 and £2 in (129]) are replaced by 

x—- = -3Ai + x£ ,i, z-r- = a:Qo,i7 
ax ax 

= -3v + x 2 Q Q , 2 , 32 

ax 

where 

~ 2^i6 2x^ 2 

<?0,1 = <?0,1 72 1" —73-^0,2 - 8XM Q ,1 + 4XW Qo,2, 

^2 ^2 

Qo,i = g (1/0,2 + ^0,2) j Qo,2 = g (2j7o,2 — ^0,2) • (33) 

The field equations are now in the form to apply Theorem (TJ giving, in a neighbourhood 
of x = 0, a three-parameter family of solutions of the form 

Ai=0(x), \ 2 = S + O(x 2 ), ^i = «i + 0(x), 

£i = 0(x), u = u 2 + O(x 2 ), v = O(x 2 ). (34) 

Restoring the original field variables, we have proven the existence of solutions of the 
field equations in a neighbourhood of the origin, satisfying the boundary conditions (|T5]) 
and analytic in So, a\, u 2 , r and £. □ 

Proposition 3 (Local existence of solutions in a neighbourhood of the event horizon) 
There exists a four-parameter family of local solutions of the field equations filMTll) 
near an event horizon at r = r^, satisfying the boundary conditions |77|) and analytic 
in r h , S h , a' h , u h , r and t. 



Proof Let ru > 0. Following [14], we take our independent variable x = r — and 
define the following dependent variables: 

p(x)=r, Al (x) = ^, A 2 (x) = 5(r), 

x 

a(r) 

ipi(x) = , ^ 2 (x) = u(r), 

x 

6W = fpl d " 6W = /fM^. (35) 
S{r) ar x ar 
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Again, these variables are the same as those in [13], apart from A2, ipi and £1. In terms 
of our new variables, the field equations fTTQHTT]) take the form 

dp 
x— = x, 
ax 

x— = -Ai + F hx + xQ hl , 
ax 

x— = xg h>2 , 
ax 

X— = ~1p! + — 

ax /cr 

X 



X 



dx Ai 



dx A X A 2 



x^T = -6 + 7\ 2 + xH h>2 , (36) 



where 



p p 3 £ 2 p 3 p 3 



2 



Ai _ _ 2f 

p pAiA 2 pA 



Gh,i = 

2A 2 £ 2 



^2 

pA 2 A 2 ' pX\ 

^2 (^ - 1) 



2 



P 2 



« 2 e 2 a 



rth,2 - — r-ry r — • (a ' J 

AiA 2 A 2 

To achieve the form ( 126]) required for application of Theorem [I] we make the further 
variable transformation 

A 2 ^i 

4>i = tpi Ai = Ai - JF h l , £ 2 = 6 - Pfc,2- (38) 

P 

The modified field equations then take the form 

dXi ~ x 
x— = -Ai + xy M , 
ax 

#1 7 , ~ 
= —0i + a\7/, 1 i, 

ax 

x^ = -£ 2 + xH h , 2 , (39) 
ax 

~ aF M £ 2 dT K1 2^i^ 2 ftF/^ 

y/1,1 = yhA ■ 



where 



dp Ai 9^2 AiA 2 <9£i 



?tfc,2 — %/i,2 

op Ai 
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Jh, 



2^1 , 2A 2 ^i iiG. 



h,2 



(40) 



p 2 \i p 3 p 2 

We may now apply Theorem [TJ which yields a four-parameter family of solutions in a 
local neighbourhood of x = 0, of the form 

p = r h + 0(x), \ 1 = 0(x), A 2 = S h + 0(x), 
= O(x), ^ 2 = + O(x), 



6 



+ 0(z) 



6 = O(x) 



(41) 



Restoring the original variables, we have proven the existence of solutions of the field 
equations in a neighbourhood of the event horizon, satisfying the boundary conditions 



( 1T7|) and analytic in r h , Sh, a' h , Uh, r and 



□ 



Comment From (110]) and the relation ([6]), we can show that for the solutions described 
by Proposition El we have 



(42) 



In order that r > corresponds to the exterior of the black hole - i.e. to the untrapped 
region - we must have fj,(r) > for r > 77,. Thus Proposition [3] gives rise (locally) to 
black hole solutions when [i'(r h ) > 0. This restricts the horizon data (£,r h , Sh,Uh,®h)- 
However, for each choice of > and £, there is an open set O of values of (Sh, uih, ot' h ) 
in M 3 that contains the horizon data (Sh, ^h, °4) = (1; 1; 0) corresponding to the trivial 
Schwarzschild-adS solution and with the property that for each (Sh, 0Jh, ot' h ) G O, 
Proposition |3] describes (locally) a black hole solution. The existence of this open set O 
is crucial for the global existence result of the following section. 

Proposition 4 (Local existence of solutions in a neighbourhood of infinity) 
There exists a five-parameter family of local solutions of the field equations ( G5HZ2P near 

^00, C\, d\, r 1 



r — y 00 ; satisfying the boundary conditions ( TJ5P and analytic in M , a 
and I. 



Proof In this case we take x = r 1 as our independent variable and define dependent 
variables as follows (again, these mostly follow [14]): 

„2~ 



Xi(x) = 2m(r) = r 



1 - M r ) + « 



X 2 (x) = S(r) 



ipx(x) = a(r), i) 2 (x) = u(r), 
2 da 2 du 



dr ' ; dr ^ 

The field equations (TTOlTTT]) . when written in terms of our new variables, take the form 



X - ^"»00,l) 



a; 



dx 
dx 



d^2 i n 
x-^r- = x y, 



-X& 



x 



dx 
dx 



00, 2, 



x- 



xH 



00,1 j 



x 



dx 
dx 



(44) 
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where 



Goo,l = "I " ^| + ^ (2 - ^) - t| - 2^00^, 
00^2 ^2 



2^| 
^A 2 



^ V / lv y 2 2 

t/00,2 — ~ -r-r, ^ ^^2^2) 



^ 2^^ 2 2 X 3 eigoo,2 

rtoo.l — ^ 1 : : 

J~ 00 <^2 



^ 2 ^ 2 (l-^2 2 ) Xe 2 (3xA!-2) S 3 &goo,l S 3 6&o,2 
oo A 2 •'00 •'00 •'00 A 2 



(45) 

No further transformations are required as the equations (1441) are already in a form 
suitable for the application of Theorem [IJ The theorem then gives us a five-parameter 
family of local solutions in a neighbourhood of x = 0, of the form 

\ 1 = M + 0(x), A 2 = l + 0(x 4 ), vj) 1 = a O0 + 0{x), 

^ 2 = co 00 + 0(x), £i = -d 1 + 0{x), £ 2 = -c 1 + 0{x). (46) 

Restoring the original variables, we therefore have proven the existence of solutions of 
the field equations in a neighbourhood of infinity, satisfying the boundary conditions 
(120]) and analytic in M, a^, Co^, c\, di, r -1 and £. □ 



4. Global existence of black holes and solitons 



In this section we turn to the issue of "globalizing" the existence and uniqueness results 
of the previous section. We consider first the black hole case and then, more briefly, the 
soliton case. 



4-1. Global existence proof for black holes 

In order to construct asymptotically adS dyonic black hole space-times, we must 
prove the global existence of solutions to the EYM equations ( TTOlTTTj) with appropriate 
boundary conditions. In this context, global existence refers to existence for all r > r^, 
where r h is the horizon radius such that fi(r h ) = 0. For each value of the adS radius 
of curvature £, the existence of a four-parameter family of such solutions is proven in 
this section using what is essentially a non-linear perturbation argument that relies on 
an application of the (Banach space) implicit function theorem - see, for example, [T5] . 
See also P3] for similar applications of this theorem to ODEs. 

Theorem 5 (Implicit function theorem) 

Let C, Y, Z be Banach spaces. Let U be a neighbourhood ofc in C, V be a neighbourhood 
°fyo in y and let F : U xV — >■ Z be a (Frechet) C 1 mapping. Suppose that F(c , y ) = 
and that the linearization d y F(c Q , yo)(0, •) : Y — > Z is an isomorphism. Then there is 
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a neighbourhood U of c in C , a neighbourhood V of y in Y and a unique continuous 
mapping G : U — >■ V such that y = G(c ) and F(c, G(c)) = for all x G U. 

Our aim then is to set up the field equations and boundary conditions in the form 
of a Banach space mapping. In this approach, C will correspond to the set (Banach 
space) of allowed initial values of the field variables, and Y will correspond to the 
space of the field variables themselves - a space of appropriately smooth functions. The 
field equations are then written in a form that corresponds to a mapping F, with the 
condition that F(c, y) = if and only if y satifies the field equations with initial values 
c. Then the trivial solution, Schwarzschild-adS, corresponds to a pair (c , yo) G C x Y 
with -^(cojyo) = 0. The existence of the mapping G of the theorem then gives rise to 
non-trivial solutions y = G(c) of the Banach space mapping, which will correspond to 
non-trivial solutions F(c,G(c)) = of the field equations and initial conditions. 

The key to applying the implicit function theorem is to write the field equations in 
a form in which the linearization about the trivial solution has certain nice properties 
throughout the interval [r h ,+oo). However, this is essentially impossible, as the 
linearized equations are singular at r = r^. To avoid this difficulty, we use the local 
existence results at the horizon (Proposition [3]) to construct solutions on an interval 
of the form [r^,?^ + 5} where 5 > 0. Then, with a suitable choice of dependent and 
independent variables, the linearized field equations do have the required properties 
throughout [r h + 5, +oo), and the implicit function theorem can be applied to produce 
solutions on this interval. This is carried out in Lemma EJ These are a continuation 
of the local solutions, and so overall we have solutions on [r^+oo) - that is, global 
solutions. The proof of this continuation is given in Proposition [7] and Lemma El 

We begin by considering again the formulation of the field equations used for the 
proof of local existence at infinity (see Proposition [4]). We recall that x = r~ l and we 
use the notation 



2 da 2 dui \ 



T 



2m 1 S 1 a 1 u 1 r 2 —,r 2 —j . (47) 



Then the equations fTT014TT]) take the form 

— = A(x, y) = (Goo.i, a; 3 £oo,2, -2/5, -y&, ^00,1,^00,2) • (48) 
The coefficients on the right hand side are given in the proof of Proposition H] (|45[) . If 

f = T QO = x 2 - Vl x 3 + i = ^ (49) 



r 



and y 2 = S are non-zero, this system of equations is equivalent to 
Q 1 :=fy 2 2 ^ + H 1 (x,y) = 0, 

Q 2 :=f 2 y2^ + H 2 (x,y) = 0, 

Q 3 :=^ + H 3 (x,y) = 0, 



Existence of Einstein- Yang- Mills dyons 15 

Q 4 :=^i + tf 4 (x,y) = 0, 
ax 

Q 5 := f 2 y 2 2 ^ + H 5 (x,y) = 0, 

Q 6 :=f 2 y 2 2 ^ + H 6 (x,y) = 0, (50) 

where 

Hi(x, y) = /y 2 2 + 2yly 2 A - fy 2 2 yj(2 - y\) + jy\ + 2/ 2 y 2 2 y 6 2 , 
H 2 (x,y) = 2x 3 (y 2 3 yl + f 2 y 2 2 y 2 6 ), 
H 3 (x,y) = y 5 , 
H 4 (x,y) = y 6 , 

H 5 (x, y) = 2fy%y 3 yj + y 5 H 2 , 

H G (x, y) = - 2/32/4 - /2/22/4(l - 2/1) - xfy 2 2 {?>xyi - 2)y 6 + 2/ 6 (x 3 i^i - H 2 ). 

(51) 

We note that the functions Hi, 1 < i < 6 are analytic for all x G R and y G M 6 . 

Let x* > — . We now introduce the Banach spaces C = M 6 , Y = C^QO, x*], M 6 ) and 
Z = R 6 x C°([0,x*],lR 6 ) and define the mapping 

F : C x Y Z : (c,y) 6CxF4 F(c,y) = f ^ ~ } C J , (52) 

where Q(x, y) = (Qi, . . . , Q§) T . Note that F(c, y) = if and only if y is a solution of 
the ODEs ( 150]) satisfying the boundary conditions y(x*) = c. We define 

c = y = (27^0,1, 0,1,0, Of. (53) 

These values correspond to the trivial solution 

H=l--^- + p, S=l, a = 0, oo = l (54) 
of the field equations (and so -F(c ,yo) = 0). The linearization d y F of F at (c ,yo) is 

v(x*) 

*Q(«,yo»v)) 



dyi f, (c o ,yo)(0,v)= I u I, (55) 



where 



SQi(x, y , v) = -^(2, y )^ + ^ ^(2, yoH, (56) 
with Vi — A straightforward calculation yields 
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bQ A = — + v 6 , 
ax 

SQs = fS-^ + 2f v 3 , 

SQ e = / 2 ^ + 2f v A + 2x(l - 3m x)f v 6l (57) 

where 

fo = x 2 - 2m x 3 + -. (58) 

Note that by the definition of x* and r^, we have fo(x) > for all x £ [0,as*]. 
The next step is to determine if the mapping 

d y F(c o ,y o )(0,-) : vey^d w F(co,yo)(0,v) eZ (59) 

is an isomorphism. This is the case if for every del 6 and every z £ C°([0, 2*], M 6 ) 
there exists a unique v £ C 1 ([0, as#], M 6 ) such that 

v(x*) = d, <5Q(x,y , v) = z. (60) 

This follows immediately by the standard theory of linear ordinary differential equations 
[T7] . The linear equations SQ = z are regular throughout the interval [0,x#] and the 
terminal conditions v(x*) = d pick out a unique solution. The fact that fo > 
throughout [0, a;*] is crucial for this argument. Thus applying the implicit function 
theorem yields the following. 

Lemma 6 There is a neighbourhood U of Cq £ M 6 and a neighbourhood V of yo £ 
C 1 ([0, x*], 1R 6 ) such that for each c £ U , there is a solution y £ V of the equations ( f50j) 
on [0,x*] toit/j y(x*) = c. The neighbourhoods U and V may be chosen so that \i > 
and S > throughout the interval [0, 2*] and i/iere is a unique y £ V for each c £ U. 

Proof We take the Banach spaces C, Y, Z and the mapping F as in the preceding 
paragraphs. Then Theorem [5] yields neighbourhoods U, V of Co,yo respectively and 
a unique continuous mapping G : U — > V with G(c ) = y and F(c, G(c)) = for 
all c £ U . By the definition of F, each y = G(c) is a C 1 solution of ( 150]) on [0, 2*] 
that satisfies the condition y(x*) = c. Next, consider an e— ball B e (yo) C V", and let 
y £ B e (y ). Then (using the standard maximum norm for C 1 ^,^*]) 

1 2/i - 2/i( c o)| < e for all as £ [0,as*], (61) 



that is 



1 2m — 2mo| < e for all r > r* > r^. (62) 

It follows that 

2m r 2 2m + e r 2 
/i(r) = 1 ^ + — >1 h tor all r >r\> r h . (63) 



Recalling that 



lM>(r h ) = 1 - ^ + i = 0, (64) 
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it is straightforward to show that a sufficiently small e > can always be chosen so that 
2m + e r 2 

1 h — r > for all r > r* > r^, (65) 

r l l 

from which it follows that /x > throughout [0, x*j. A similar, but more straightforward 
argument shows the same conclusion for S. We now take V = B e (yo) for the sufficiently 
small positive e found above, and take U = G~ l {V). By continuity, this is open (and 
must contain c ). Positivity of \x and S show that / = x 2 — x 3 U\ + 1~ 2 and y 2 = S are 
bounded away from zero on the interval [0, x*]. Then uniqueness of y £ V as a solution 
of flSTJj) for a given c £ U follows by the standard argument of estimating the growth 
of 8y, the difference of two potentially distinct solutions. (It is crucial to have bounds 
on the coefficients of the derivatives in (1501) for this standard argument to apply; these 
coefficients are powers of / and 2/2 ■) D 

Comment Lemma[6]establishes the existence of solutions of (loTJj) on the interval [0, x*]. 
These solutions have / = /i/r 2 >0 and y 2 — S > 0, and so correspond to solutions of 
(JTOHTTj) , defined on an interval of the form [r*, +00) where r* > r^. In the following 
results, we prove that an open subset of these solutions provide a continuation of the 
local solutions at the horizon, and thus yield global solutions defined for all r > r^. This 
proof relies on the geometric structure of the set of local solutions near the horizon. 

Proposition 7 (Global existence of black hole solutions) 

Let rh > 0,A < and define i = — ^. Then there is an open neighbourhood T hh of 
go = (1,1,0) £ 1R 3 such that for each g = (S^, uj h , a' h ) £ r bh , there exists a global 
solution of the EYM field equations 001471]) defined on [r h , +00) with the asymptotic 
behaviour [T7\ ) as r — > r£ and the asymptotic behaviour (E2J) as r — > 00. These solutions 
are differentiate on [rh, 00) (C 1 for fi, S and C 2 for a,u) and satisfy /i(r) > for all 
r > Th- Since fi(rh) = 0, the region r > r^ corresponds to the exterior of a black hole. 

Proof By Proposition [3] and the comment which follows its proof, there is a 
neighbourhood T of g £ M 3 such that for each g £ T, we obtain a solution of the 
EYM field equations (TlOHTT]) on an interval [r^, r^ + S) corresponding locally to a black 
hole solution. These solutions are analytic in r and in g on this interval. Now choose 
x* = r~ l with r/j < r* < + 5. Then it follows by the same proposition that the 
solution vector y defined by ( 1471) is analytic in g at x = x*. Let y(x; g) denote the 
solution y(x) with horizon data g and define a mapping 

(P-.T^R 6 : g£r^0(g) = y (x,;g). (66) 

Note that <fi is analytic on V as a consequence of Proposition [3j The image M. of 0, 

M = 0(r) = {0(g) £ R 6 : g £ T} (67) 

is a 3-parameter subset of M 6 that contains the point Co = (2mo, 1, 0, 1, 0, 0). (Recall 
that this vector c £ M 6 corresponds to the trivial solution of the field equations with 
horizon data go-) We can show that without loss of generality, A4. is a three-dimensional 
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submanifold of M 6 in an neighbourhood of c : this is proven in Lemma M below. Now 
define U = Ai n U, where U is the neighbourhood of c described by Lemma [6j By 
this lemma, for each c £ U, there is a unique solution of the EYM equations on [r*, oo). 
Since Ai is a manifold, it follows that U is an open subset of Ai. Thus each g £ r bh 
gives rise to a solution of the EYM equations on [rh, +00), where 

T bh = <p-\U) = {g £ r : y(x*; g) £ U} (68) 

is an open subset of T, being the inverse image of an open set with respect to an analytic 
(and hence continuous) mapping. The asymptotics carry over from the corresponding 
local results at the horizon and at infinity respectively. □ 

Lemma 8 Let T and be as in the proof of Proposition^ Then there is an open subset 
Vi ofT containing g such that Ai\ = (j>(Ti) is a three-dimensional submanifold 0/IR 6 . 



Proof Let us write g = (£, r\, () = (Sh, <^h, Oi' h ) F for the horizon data, and write 

M = {^, f/l ()ei 6 :(e,^()Gr}. (69) 

We show that in a neighbourhood of 0(go), we can describe Ai as the graph of a C 1 
function. To see this, we recall that the local solution described by Proposition |3] is 
analytic in both r and the horizon data g for values of r sufficiently close to r^. The 
quantity x* corresponds to a value in this range, and so <fi is analytic in g. Recalling 
the definition (147j) and the asymptotics derived as part of Proposition [31 we have 

02 = e + 0(r-r h ), 
04 = 7) + 0(r - r h ), 



05 = r 2 h C + 0(r - r h ). 



(70) 



In fact, by analyticity, the remainder terms can be replaced by convergent power series 
in r — rh, and so may be differentiated with respect to £, r] and £. That is, we can apply 
the rule 

d 



[0(r - rh)] = 0(r - r h ), 



and the corresponding results for rj and (. The coefficients of these power series are 
analytic functions of (£, 77, C)- Thus the Jacobian matrix of the mapping : (£, 77, C) — > 
(02, 04, 05) has the form 

/ l + 0(r-r h ) 0{r-r h ) 
0(r-r h ) l + 0(r-r h ) 
0(r-r h ) 0(r-r h ) 



(9(02, 04, 5 ) 



with determinant 
J(0) 



V 



0(r - r h ) 
0(r - r h ) 
r 2 h + 0(r-r h ) J 



,(71) 



<9(02,04,05) 



r\ + 0(r - r h ). 



(72) 
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Thus for chosen sufficiently close to r^, this determinant is non-zero and the inverse 
function theorem [15] shows that there exist C l functions V'1,'02,'03 and an open 
neighbourhood I\ of (1, 1, 0) such that 

\ 



(yA 




( 


MZ>v,0 ^ 




(A 


( 


2/4 










V 




\y 5 J 






h(£,v,Q J 






\ 



^112/2,2/4,2/5 J 

^2(2/2,2/4,2/5) 



(73) 



M l 



■ (2/2,2/4,2/5) e 0(ri 



(74) 



for (£, 77, C) G Ti and (2/2,2/4,2/5) £ ^Fi). Now define .Mi = <f>(Fi). This can be written 

as 

/ 01(2/2,2/4,2/5) \ 
2/2 

^3(2/2,2/4,2/5) 

2/4 
2/5 

V 06(2/2,2/4,2/5) / 

where ^ = ^0^1- that is, ^(2/2,2/4,2/5) = 0i(^i (2/2, 2/4, 2/5), ^1(2/2, 2/4, 2/5), ^1(2/2, 2/4, 2/5)), 
z = 1, 3, 6. This shows that .Mi is the graph of a smooth (C 1 ) function, and is therefore 
(counting dimensions) a three-dimensional submanifold of M 6 - again, by an application 
of the inverse function theorem. See, for example, chapter 3 of [18] for further details 
of this last point. □ 



(Vi \ 




2/2 




2/3 




2/4 




2/5 





Comment Proposition [7] has proven the existence of non-trivial black hole solutions of 
the field equations (ITDHTTT) in a neighbourhood of the trivial Schwarzschild-adS solution. 
The proof is valid for any value of the adS radius of curvature £, but does not give us any 
information about the size of the neighbourhood in which we have non-trivial solutions. 
Numerical investigations j6j [10] indicate that the size of this neighbourhood will shrink 
as £ increases (or, as the magnitude of the cosmological constant |A| decreases). Since 
the results of section [3] tell us that the solutions are analytic in the horizon data, non- 
trivial solutions which are sufficiently close to the trivial solution will be such that the 
magnetic gauge field function u(r) has no zeros since the trivial solution u(r) = 1 has 
no zeros. We anticipate that at least some of these solutions will be stable. 



4-2. Global existence proof for solitons 

The method applied above to establish the global existence of asymptotically adS black 
holes carries over with minimal changes to the soliton case. We highlight those changes 
here and then quote the relevant result. 

The argument proceeds unchanged up to and including the definition ( 1521) of the 
Banach space mapping F . The linearization is now about the trivial data point and 
trivial solution 

co = y = (0,1, 0,1,0, Of (75) 
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corresponding to pure anti-de Sitter space-time wherein 



/i(r) = l + -, S{r) = l, a(r) = 0, u(r) = 1. 



(76) 

KT 

The linearization fl55|) maintains the same form fl57|) but with mo = 0, and hence 
fo = x 2 + £~ 2 . Lemma [6] thus carries through. In the present the reciprocal 

of a value r* with < r* < 5, where the local existence result Proposition [2] applies on 
[0,5). These local solutions satisfy the boundary conditions ( IT5l) . which are equivalent 
to 

S + S 2 r 2 + O(r 3 ), 
= l + w 2 r 2 + 0(r 3 ), 
2w 2 r 3 + 0(r 4 ). 



2/i 

2/3 
2/5 



2m 



r 2 a' 



= 2m 3 r 3 + 0(r 4 ), y 2 
a x r + a 3 r 3 + 0(r 4 ), 
= a x r 2 + 0(r 3 ; 



S 



2/4 



2/6 



rV 



(77) 



The free parameters here are (£,?7,C) = (So, «i, ^2), with the other parameters being 
given in terms of these (see ( 115}) ). These represent the free initial data at the origin for 
solitons, and local existence of an analytic solution applies for data in a neighbourhood 
T of the trivial data (£,?7,C) — (1)0,0). As in the black hole case, the local existence 
theorem gives rise to an analytic mapping defined as in fl66|) . Corresponding to 0, we 
define the mapping : (£, r], () (0 2 , 03, 

/ l + 0(r 2 : 



J 4), 



which has Jacobian matrix of the form 



•9(02,03,04) 



with determinant 
J(0) 



V 



0(r 3 ) 
0(r 3 ) 



0(r 2 ) 
r + 0(r 3 ) 
0(r 3 ) 



0(r 2 ) \ 
0(r 3 ) 
r 2 + 0(r 3 ) y 



^(02,03,04 



r 3 + 0(r 6 ). 



(7? 



(79) 



^,^/,C) 

Thus for sufficiently small r, we see that is a diffeomorphism in a neighbourhood of 
the trivial data. This is the last technical result required to ensure that the analogues 
of Proposition [7| and Lemma [8] carry through to the soliton case, and so we can quote 
the following result. 

Proposition 9 (Global existence of soliton solutions) 

Let A < and define = — -g. Then there is an open neighbourhood T so1 of 
go = (1, 0, 0) e R 3 such that for each g = (S , a±, uj 2 ) G r so1 , there exists a global solution 
of the EYM field equations (T(MTl\) defined on [0, +oo) with the asymptotic behaviour 
/ fl5j) as r — > and the asymptotic behaviour [2D\) as r — > +oo. These solutions are 
differentiable (C 1 for fi and S , C 2 for a, go) on [0, oo) and satisfy fi(r) > for all r > 0. 



Comment Proposition |9] has proven the existence of non-trivial soliton solutions of 
the field equations ffTOTfTTT) in a neighbourhood of the trivial pure adS solution. As with 
the corresponding result for black holes, the proof is valid for any value of the adS radius 
of curvature £, but does not give us any information about the size of the neighbourhood 
in which we have non-trivial solutions. For solitons as well as black holes, numerical 
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work [BJ H3] shows that the size of this neighbourhood will shrink as I increases. Since 
the results of section [3] tell us that the solutions are analytic in the initial data at the 
origin, non-trivial solutions which are sufficiently close to the trivial solution will be 
such that the magnetic gauge field function u)(r) has no zeros, solutions which are of 
particular interest because we expect at least some of them to be stable. 

5. Conclusions 

In this paper we have studied dyonic soliton and black hole solutions of the su(2) EYM 
equations in adS. These solutions have been found numerically some years ago [6] but, 
while the purely magnetic system has been studied in depth analytically, the dyonic 
system has received comparatively little attention in the literature. Here we have proven 
analytically the existence of solutions of the field equations representing dyonic solitons 
and black holes. The approach is two-fold: firstly, we prove the local existence of 
solutions near the origin, black hole event horizon and infinity, where the governing 
differential equations are singular. Global existence of solutions is then proven using a 
non-linear perturbation argument which relies on an application of the Banach space 
implicit function theorem. This approach to proving global existence of solutions of the 
EYM equations in adS is novel - previous work (for example [5]) has used a "shooting" - 
type approach based on showing that the local solutions in a neighbourhood of the origin 
or event horizon can be extended to large values of r, and then matching onto the local 
solutions in a neighbourhood of infinity. 

The result proved in this paper is valid for any non-zero adS radius of curvature i 
(or, equivalently, any negative cosmological constant A < 0). Since it is known [11] that 
there are no non-trivial dyonic solutions of the EYM field equations when the space-time 
is asymptotically flat rather than asymptotically adS, a natural question is how our proof 
breaks down in the A = 0, i — > oo case. The key difference between the asymptotically 
flat and asymptotically adS cases is the boundary conditions as r — > oo. In particular, 
Proposition |4] (giving local existence of solutions near infinity) does not hold as currently 
stated because Too = when x = and £ — > oo, so the ordinary differential equations 
( T44|) are not regular. This result breaks down because, in asymptotically flat space, the 
boundary conditions at infinity imply that either oorw must vanish as r — > oo. 

The global existence result in this paper proves the existence of non-trivial solutions 
in a neighbourhood of the trivial solution a(r) = 0, u(r) = 1, M(r) = constant, 
S(r) = 1, corresponding to either pure adS or the Schwarzschild-adS black hole. In 
particular, at least some of these solutions, sufficiently close to the trivial solution, are 
such that the magnetic gauge function u(r) has no zeros. These solutions are particularly 
interesting because one may conjecture, in analogy with the purely magnetic case [5j E] 
that at least some such solutions will be linearly stable. We will return to this question 
in a forthcoming publication. 

It might also be of interest to investigate whether the analytic approach taken in 
this paper can be applied to soliton and black hole solutions of Einstein gravity coupled 
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to other matter fields. For example, adding scalar fields to the model (possibly with 
coupling to the gauge field) would be relevant to supergravity models with non-abelian 
gauge fields [19]. We leave such questions for future research. 
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